Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 22 October 2009 



(MN style file v2.2) 



Evolution of supermassive stars as a pathway to black hole 
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ABSTRACT 

Supermassive stars, with masses > 10^ Mq, are possible progenitors of supermassive 
black holes in galactic nuclei. Because of their short nuclear burning timescales, such 
objects can be formed only when matter is able to accumulate at a rate exceeding 
~ IMq yr~^. Here we revisit the structure and evolution of rotationally-stabilized 
supermassive stars, taking into account their continuous accumulation of mass and 
their thermal relaxation. We show that the outer layers of supermassive stars are not 
thermally relaxed during much of the star's main sequence lifetime. As a result, they do 
not resemble n — 3 polytropes, as assumed in previous literature, but rather consist of 
convective (polytropic) cores surrounded by convectively stable envelopes that contain 
most of the mass. We compute the structures of these envelopes, in which the equation 
of state obeys P/p^/^ oc ]Vp/'^{R), where M{R) is the mass enclosed within radius R. 
By matching the envelope solutions to convective cores, we calculate the core mass 
as a function of time. We estimate the initial black hole masses formed as a result of 
core-collapse , and their subsequent growth via accretion from the bloated envelopes 
( "quasistars" ) that result. The seed black holes formed in this way could have typical 
masses in the range ^ lO"* — IO^Mq, considerably larger than the remnants thought to 
be left by the demise of Pop III stars. Supermassive black holes therefore could have 
been seeded during an epoch of rapid infall considerably later than the era of Pop III 
star formation. 



Key v^rords: black hole physics 
quasars: general 
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1 INTRODUCTION 

Supermassive stars were proposed by Hoyle & Fowler 
(1963a,b) as a means to meet the prodigious energy re- 
quirements of radio galaxies (Burbidge 1958), and slightly 
later as a model for quasars. Hoyle & Fowler recognized 
that such objects could not persist for much more than a 
million years, if sustained by hydrogen burning, and sug- 
gested that they might collapse to form black holes once 
their nuclear fuel was exhausted. The pulsational stability 
of supermassive stars is a crucial issue, because they are ra- 
diation pressure-dominated and therefore have an adiabatic 
index very close to 4/3 (which yields neutral stability to ra- 
dial pulsations for a Newtonian, self-gravitating body with 
no rotation). Small general relativistic corrections have a 
destabilizing effect, preventing nonrotating stars more mas- 
sive than a few xlO^ Mq from attaining a phase of sta- 
ble hydrogen burning before collapsing (Iben 1963; Fowler 
1964). However, a dynamically insignificant level of rotation 
— especially differential rotation — can stabilize stars as 



massive as ~ or more (Fowler 1966). Additional sta- 

bilizing effects due to magnetic fields and turbulence were 
later considered by other authors (e.g., Bisnovatyi-Kogan, 
Zel'dovich & Novikov 1967; Ozernoy & Usov 1971). 

In this paper, we are interested in whether stable super- 
massive stars might be precursors of seed black holes that 
eventually grow to large enough mass to power quasars and 
populate the nuclei of present-day massive galaxies. As we 
will argue in § 5, the holes created inside a supermassive 
star are likely to have masses of a few percent of the star's 
final mass. Thus, to obtain seed black holes with masses 
~ 10* — 10^ M0, comfortably larger than the seeds probably 
left behind by the collapse of Population HI stars, we need 
to consider supermassive stars with masses > IQ^'Mq. We 
will focus on this mass range in our analysis □ 



^ Note, however, that the term supermassive star is often used 
to refer to the mass range M. > 5 X lO^M© (Fuller, Woosley & 
Weaver 1986). 
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Once the existence of a stable supermassive star is pos- 
tulated, its structure would appear to be quite simple. In or- 
der for radiation pressure to support the star against gravity, 
the radiative flux seen by every element of mass in the star 
must equal the local Eddington limit, Fe = GM{R)c/ nr'^ , 
where k is the local opacity and M(R) is the mass enclosed 
within R. The opacity inside a supermassive star is likely to 
be roughly uniform, with electron scattering dominant ev- 
erywhere, but M{R) is a monotonic function of radius. Since 
the entire luminosity is produced by thermonuclear reactions 
within a small region of the core, the star must be highly 
convective with the fraction of the total luminosity carried 
by convection varying as Lcom- / itot = 1 — M (i?) /A/, , where 
Mt is the total mass of the star. If the convection is cfBcicnt, 
then a supermassive star should bo described accurately by 
an n = 3 (7 = 4/3) polytrope, according to the Lane-Emden 
equation. 

Such models, however, beg the question of how a su- 
permassive star might realistically form. Creating a 10® 
star requires the very rapid accumulation of gas. Since the 
thermonuclear timescale for a star burning hydrogen at the 
Eddington limit is ~ 2 Myr, independent of mass, and will 
be less if only part of the fuel is burned, the infall rate 
required to create a supermassive star of mass must 
exceed 0.5{M^/1Q^Mq) Mq yr"^ This is orders of mag- 
nitude larger than the rate at which matter condenses to 
form normal stars within molecular clouds, or even the rate 
(~ 1O~^M0 yr~^) at which matter came together to form 
Pop III stars in pregalactic dark matter halos. 

The maximum rate at which matter can collect is given 
by ~ v^/G = 0.2(t;/10 km s"^)-^Af0 yr"\ under the as- 
sumption that the gas is self-gravitating and collapses at 
its free-fall speed v. In effect, this means that supermas- 
sive stars can only form in systems with virial tempera- 
tures exceeding 10* K. There is considerable controversy 
over whether such high rates of inflow can occur without 
most of the gas fragmenting and forming stars before reach- 
ing the center. Despite early suggestions that fragmentation 
is avoidable only if the gas temperature remains close to the 
virial temperature (Bromm & Loeb 2003; Begelman, Volon- 
teri & Rccs 2006) — and thus only in systems lacking both 
metals and molecular hydrogen — recent simulations sug- 
gest that fragmentation may be suppressed even when the 
gas is much colder (Wise, Turk & Abel 2008; Levine et al. 
2008; Regan & Haehnelt 2009), possibly due to the con- 
tinuous generation of supersonic turbulence (Begelman & 
Shlosman 2009). 

Conditions of rapid infall are ideal for producing the 
high levels of entropy required inside supermassive stars. 
Indeed, simple estimates show that the initial entropy of 
the gas joining the star cam be much larger than the 
value required for equilibrium. We argue that the nature 
of radiation-pressure support provides a mechanism for au- 
tomatically regulating the entropy of the gas joining the su- 
permassive star, so that it is neither too large nor too small. 
This is due to the existence of a "trapping radius" within the 
infalling gas, outside of which radiative diffusion can release 
excess entropy. Gas with too much entropy will be forced to 
expand until it gives up the excess. However, this does not 
imply that supermassive stars should have uniform speciflc 
entropy and therefore resemble polytropes. 

In this paper we study the evolution and fate of su- 



permassive under the assumption that they grow by rapid, 
continuous infall. In §§ 2 and 3, we argue that gas joins the 
star with increasing specific entropy as a function of time. 
Since the Kelvin-Helmholtz timescale is longer than the age 
of the star for early times and high infall rates, much of 
this this entropy stratification is preserved over the life of 
the star. We therefore conclude that supermassive stars are 
not necessarily well-represented by n = 3 polytropes, but 
rather can have a more complex structure with a convec- 
tive (polytropic) core surrounded by a convectively stable 
envelope that contains most of the mass. The entropy in 
the envelope roughly satisfies P/p*^^ oc M^^^{r) (we term 
this entropy law "hylotropic" ) . Hydrogen burning in the core 
starts when the star's mass and entropy are both relatively 
low, and adjusts to sustain the star through its more massive 
stages (§ 4). 

In § 5, we discuss the formation of the seed black hole 
after the exhaustion of core hydrogen, and its subsequent 
growth inside the remnant of the supermassive star. Because 
of the likely importance of rotation in the collapsing core, 
only a small fraction of the supermassive star collapses to 
a black hole initially. The energy liberated during the for- 
mation of the black hole inflates the remainder of the star 
into a bloated object that resembles a red giant, which we 
have previously termed a "quasistar" (Begelman et al. 2006; 
Begelman, Rossi & Armitage 2008). The black hole contin- 
ues to grow by accretion from the quasistar envelope, until 
the photosphcric temperature drops to the point where the 
quasistar undergoes an "opacity crisis" and disperses under 
the influence of radiation pressure. We summarize our re- 
sults in § 6 and comment on the possible contribution of 
supermassive stars to the ionizing radiation field at high 
redshifts. 



2 SELF-CONSISTENCY OF FULLY 
CONVECTIVE MODELS 

A thermally relaxed supermassive star must be fully convec- 
tive and therefore well-represented by an n = 3 (7 = 4/3) 
polytrope (Hoyle & Fowler 1963a). The structure is deter- 
mined by the solution, 63{x), of the Lane-Emden equation: 

with boundary conditions ^3(0) = 1, ^3(0) = 0, where 
X = R/R:i is the radius normalized to Ri — (Pcji^Gfi^^^'^ . 
Pc and pc arc the central pressure and density, respectively, 
with the pressure and density elsewhere being given by 
P{x) = Pc03{x) and p{x) = pc93{x). Density is related to 
pressure via an entropy parameter K = Pj p^^^ , which de- 
pends uniquely on the mass of the star: K = 0.364GMr^ 
The dimensionless radius of the star is given by x» = 6.897, 
corresponding to 

i^, = 5.8 X lO^mi^e Tjyg cm, (2) 

where m.,6 = M,/10^Mq and Tc = lO^T^.s K is the cen- 
tral temperature. Rather than assuming that the star sim- 
ply exists with a fixed mass, we assume that the star's 
mass is growing at a rate M, = m* Mq yr~^, so that 
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Given the growth rate of the star, we can assess whether 
the star is thermaUy relaxed by comparing 7?« to the "trap- 
ping radius," Rn, which separates the outer region where the 
infaUing gas leaks radiation, from the inner region where the 
flow is approximately adiabatic (Begelman f978). For elec- 
tron scattering opacity with k — 0.34 cm^/g, we have 



mass, P/p' 



4/3 



Af2/^ Since mass increases with radius, this 



Rtv = -^^^ = 5.7 X fO^^rh, cm. 
47rc 



(3) 



Suppose R* lies inside the trapping radius. In this case, the 
infaUing gas must already be adiabatic by the time it joins 
the star. This is problematic, because the specific entropy of 
the gas in an n = 3 polytrope (as measured, e.g., by K) is 
uniquely determined by the mass. Matter about to join the 
star will generally have much higher entropy than the re- 
quired value. (This can be ascertained, in a specific case, by 
comparing the ram pressure with the density of gas in free 
fall — but it would be true in most circumstances.) There- 
fore, infaUing matter must lose entropy in order to join a 
polytropic star. But if R, < Rti, the gas behaves adiabat- 
ically near the star and such entropy loss is impossible. In 
this case the star cannot be thermally relaxed. On the other 
hand, thermal relaxation is possible if J?* > i?tr. 

Comparing equations Q and we see that super- 
massive stars of a given mass can become thermally relaxed 
only after a certain amount of time has passed. 



(4) 



or equivalently, for a given accretion rate, if tMyr > m^T^^g. 
The "main sequence" lifetime for a fully convective star is 
« 2 Myr, where we have assumed a hydrogen mass fraction 
X = 0.75. This implies that no star more massive than ~ 
4 X 10^T~g Mq can be thermally relaxed during its lifetime. 



3 PARTIALLY CONVECTIVE 
SUPERMASSIVE STARS 

If matter falling onto a growing supermassive star has the 
"wrong" entropy to maintain a polytrope, what kind of 
structure results? It appears that the acquisition of entropy 
can be a self-regulating process. If gas with too much entropy 
per unit mass is added to an n = 3 polytrope, for example, 
the star will expand to the point where its radius exceeds the 
trapping radius and it can radiate away the excess entropy. 
If, conversely, gas with too little entropy is added to the 
polytrope, the star will shrink until thermonuclear reactions 
increase the entropy in the interior. Since a supermassive 
star is built up by adding gas with "too much" entropy, this 
implies that the radius of a supermassive star is of order the 
trapping radius associated with the infall rate. For constant 
M, the radius of a supermassive star is independent of mass. 
This is impossible if the entropy is uniform throughout the 
star (that would be a polytrope, with Rt oc M^^), but it is 
possible if the specific entropy increases with radius. 

At each stage in the star's growth, the entropy of the 
newly-added gas adjusts to the "correct" value for the star's 
current mass. As more matter is added, the earlier layers 
are squeezed to smaller radii, and their entropies are frozen 
in. From the homology scalings P ~ GM^/R* and M ~ 
pR^ , which must be satisfied for each layer of mass, we see 
that the entropy must be an increasing function of enclosed 



implies that the envelopes of unrelaxed stars are stable to 
convection. 

We now construct quantitative models for these unre- 
laxed envelopes, and show how they can be matched to con- 
vective cores. Because the equation of state depends explic- 
itly on the enclosed mass, it is convenient to use Lagrangian 
variables. Assuming hydrostatic equilibrium, the equations 



'i-npR^R! = 1 



P' 



GM 



(mass conservation) (5) 
(hydrostatic equilibrium) , (6) 



47ri?-* 

where a prime denotes differentiation with respect to M, to 
which we add the equation of state 

P = Ap*/^M''/^ (7) 
where j4 is a constant. 

3.1 Hylotropes 

Equation (O is a special case of a more general class of 
equations of state of the form 



P = Ap^/^Af". 



(8) 



Because of the dependence on mass, we propose to call such 
equations of state hylotropic and the derived structures hy- 
lotropes (from Gk. hyle, "matter" + tropos, "turn")0 By 
eliminating p through equations ((5]), ([6]), and ([8]), we obtain 



EL 
p 



a 
M 



8R!_ 
3 R 



4^ 
3 R' 



(47r)^/3GM^- 
A 



p/\ 4/3 



Eliminating R in favor of lu = MR' /R, we obtain 



,,ui , 3q 
M— = 1+ — 
w 4 



- 3k; + jyAf 3 



4/3 



(10) 



where r) = 3(47r)^''^G/4A. Equation fTO]) reduces to a La- 
grangian equivalent of the Lane-Emden equation for n = 3 
polytropes when a = Q. For the case of interest here, 
Q = 2/3, we have 



w 2 



(11) 



which is unique in having no preferred mass scale {rj is 
dimensionless) . Equation (llip admits power-law solutions, 
R (X with w — const., provided that the right-hand-side 
has a real root. This occurs for rj < rjait = 9/2^^^ ~ 1.7858, 
i.e., for sufficiently large entropy (since rj scales inversely 
with A). The power-law solutions are unbounded and there- 
fore unphysical, just as n = 3 polytropes are if the entropy 
is too large for the mass, i.e., if if > OMAGM,'-". 
From eq. we see that 



Ip 
w — , 

3p ' 



(12) 



where p is the mean density inside R. Since p must be a 
monotonically decreasing function of radius, we must have 
w ^ 1/3 everywhere. All solutions satisfying this constraint 

^ We thank A. Accardi and G. Lodato for suggesting this com- 
pound word and providing its etymology. 
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at some inner boundary Mm tend to power-laws at M 3> 
Mm when ry < r/crit; therefore, we are interested only in 
models with r) > rjait, all of which are bounded. 

The bounded solutions all have w increasing monotoni- 
cally with M and R; the outer edge of the envelope is reached 
where w diverges. In practice, it is easiest to treat w as the 
independent variable, computing M(w) and R{w) by inte- 
grating 



din M 

dw 

dlnR 
dw 



— (1 — 2w) -f rjw 



4/3 



(13) 
(14) 



The normalizations of M and R are arbitrary; one need only 
specify an initial value of w and integrate outward. 



3.2 Matching to a Convective Core 

Suppose that an a = 2/3 hylotropic envelope matches to 
a convective core at Mm, Rm- The equation of hydrostatic 
equilibrium implies that P and P' are continuous, while the 
continuity of the specific entropy implies that p is contin- 
uous as well. The mass conservation equation then implies 
that w — MR' /R is continuous. Only p' changes discon- 
tinuously across the boundary. Matching the hylotropic and 
polytropic equations of state at the boundary gives a rela- 
tion between the polytropic entropy parameter K and the 
hylotropic entropy parameter 77: 



K = AM^'' = 



3(47r)i/^GM, 
4r7 



2/3 



(15) 



The mass enclosed within the dimensionless matching 
radius x™ is 



Mm = 



(16) 



Substituting for K from eq. (|15|l . we find that Mm cancels 
out and we are left with a relationship between Xm and r): 











ir 


O^x^dx 




Jo 





2/3 



(17) 



From eq. (|12p . we see that the matching radius also deter- 
mines the value of w at the match point. 



1^3 (^Xn 



(18) 



Thus, by selecting a matching radius Xm we specify both 
the entropy scaling parameter of the envelope (i.e., 77) and 
the initial value of Wm to use in integrating equations H13p 
and (|14[) . The critical entropy parameter 77crit corresponds 
to a critical matching radius 1.2957 < Xcrit < 1.2958 and a 
critical value of Wm ~ Wait ~ 0.455. There are no bounded 
envelope solutions for Xm < Xait (wm < Wait), since these 
correspond to 77 < 77crit. 

As the critical values of these parameters are ap- 
proached from above, the ratios of the envelope mass to 
the core mass [M^/Mm) and the envelope radius to the 
core radius {Rt/Rm) diverge. For example, for Xm = 
1.2958, numerical integration of equations (fT3)) and (fT4)) give 
M-t/Mm = 2.1 X 10* and R,/Rm = 3.6 x 10^^, respectively. 




Figure 1. Ratios of envelope-to-core mass (lower curve) and ra- 
dius (upper curve) for hylotropic envelopes matched to polytropic 
(n = 3) cores, as a function of the dimensionless matching radius 
^m • As Xm increases from 1.3 to 2, the mass ratio decreases from 
10^ to just 2.27. Over the same interval, the entropy parameter 
rj goes from 1.80 to 3.09. 




log. 



Figure 2. Density vs. radius for hylotropic envelopes surrounding 
convective cores, normalized to conditions at the matching radius, 
for three ratios of envelope mass to core mass. 



However, these ratios decline extremely rapidly as the pa- 
rameters increase only slightly above their critical values. 
Fig. 1 shows that all interesting ratios of envelope-to-core 
mass and radius are contained within a narrow range of 
matching radii and entropy parameters. 

Envelope profiles are easily computed using equations 
(O and ((7|. Figures 2 and 3 show plots of the normalized 
density and pressure. 



M 



(Pm) 



Wm 
W 



P 
P 



\Mm) \Rm) \ W J 



4/3 



(19) 



for M^/Mm = 10,100,10^. The density decreases with ra- 
dius roughly as R~"'^^ while the pressure decreases as R"^ . 
However, a power-law is a poor approximation: the slopes 
vary continuously. This curvature on the log-log plot is more 
evident if one plots density or pressure vs. mass. 
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Figure 3. Same as Fig. 2, but for envelope pressure. 
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Figure 4. For envelope-to-core mass ratios in the range 10— 10'^, 
the ratio of envelope-to-core radius is very close to the square of 
the mass ratio. 



It is also apparent from Figures 2 and 3 that there is a 
simple approximate scaling relation between the envelope- 
to-core ratios of radius and mass; 



f— V 



(20) 



Rm 

As Fig. 4 shows, this result is not exact, but for mass ra- 
tios > 10 it is remarkably accurate. We will see below that 
this relationship allows us to deduce a simple evolutionary 
sequence for supermassive stars as they grow in mass. 

3.3 Core— Envelope Co-evolution 

As in a normal star, we expect the extreme temperature- 
sensitivity of thermonuclear reactions to regulate the cen- 
tral temperature of the core, and therefore Pc — aT^ /i, to 
a nearly constant value as the core evolves. Given that Xm 
(and therefore rf) is roughly constant as well, the relation- 
ship between mass and radius of the core becomes a simple 
scaling law: 



Rn 



3/4 



4N 1/4 



G 



l/4p-l/4,,l/2 



AC oc M, 



1/2 



(21) 



which allows us to eliminate Mm in favor of Rm in eq. (|20|l . 
Since Af, = Kht is given as a function of time, eq. (|20|) is 



readily reduced to a relation between Rt and Rm- A second, 
independent relation linking these two radii would therefore 
suffice to determine the co-evolution of the core and enve- 
lope. This relationship is provided by the virial theorem. 

According to the virial theorem, the total energy of 
a radiation-dominated star is very much smaller than the 
canonical value of — GAf^/i?*, Etot ~ —Eg — E^ot + -Egr, 
where Eg is the thermal gas energy, E^ot is the rotational 
kinetic energy, and Egk is the energy correction associated 
with general relativistic effects (Fowler 1964). The relativis- 
tic correction must be outweighed by the gas pressure and 
rotational terms in order to ensure stability, and we hence- 
forth neglect it. If we assume that fresh matter joining the 
star is "cold," so that 4P/p <^ GM-t/R*, then accretion de- 
creases the energy of the star at a rate —GMtMt/R*- In ad- 
dition, the star leaks energy at the Eddington limit, so that 
the total rate of energy loss is —GM^Ph/R* — Le{M,) = 
-{GM,Mt/R*){l + Rt/Rti), where we have used the defi- 
nition of the trapping radius from eq. The corrections 
to this energy loss due to changes in the gas or rotational 
energy are negligible (assuming that the rotation is well be- 
low Keplerian). To maintain equilibrium, nuclear reactions 
in the core must compensate for virtually all of this energy 
loss. 

The virial theorem also tells us the energy of the en- 
velope and core separately, in terms of Rm and Pm, the 
pressure at the matching radius. The envelope has a nega- 

E,ot(> Rm)-4nPmRl. 



tive energy equal to —Eg{> Rm) - -Erot(> Rm) 
where, under our assumptions, Pm ~ Pcdiixm) is approxi- 
mately constant. [The core energy is —Eg{< Rm) — £rot(< 
Rm) + iirPmRm-] The envelope loses energy at Rt and 
gains energy at Rm by a combination of radiative diffusion, 
at a rate comparable to the Eddington limit for the core 
mass. Mm, and advection of energy across the core-envelope 
boundary. Since 

PmRm is of order GM^/R,, we can neglect 
the change of gas and rotational energy, and we also crudely 
estimate the energy per unit mass crossing the boundary to 
be GMf/Rt (i.e., the same as at i?,). Matching the change 
of envelope energy to the energy loss at the surface, we ob- 
tain 

An^iPmRt) = 
at 



GM,M, 
R, 



Mm 

A/. 



+ 



R, 
Ru 



Mm 

M^ 



(22) 



If we consider Mm/M, to be a slowly varying function of 
time (we will see below that it scales roughly as t^^^), we can 
further approximate Mm/M, ~ Mm/M^, which gives us the 
desired behavior that the left-hand side vanishes when Mm 
approaches Af», and also should be accurate when Mm <S 
Af.. 

Using equations (I21|l and (|20p . and assuming constant 



(23) 



Mt, Pc, and Xm, we then obtain the approximate equation 

t dtp* \ M, ) \R, RtJ ' 
where x = ("27 /4rj^)Xm63{xm) ■ We then obtain 
2x+^~l kM. 



R 



2X + 



Af. 

Mm 



-Rti 



Mm 



(24) 
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As we guessed from our simple analysis based on radiation 
trapping arguments, the radius of a partially convective su- 
permassive star is roughly proportional to the accretion rate. 
Equation (|24p . and eq. (|23p from which it is derived, are 
only approximate because xi^m) does vary with Xm, which 
in turn varies weakly with M-t/Mm (Fig- !)• But these re- 
lations show that R, ^ (0.4 - 0.8)i?tr for Mm/M, > 0.1, 
and is no larger than « 1.2i?tr for Mm/M, > 10"^. For 
the purpose of the estimates to follow, it will suffice to take 
R*/Rti ~ 0.6 for a partially convective star. 

Once we have an expression for it is straightfor- 
ward to calculate the evolution of the core. In terms of the 
normalized quantities defined earlier, the core mass is given 
by 

2/3 



■ 6.2 X 10' 



7,1/2 (7?,/i?t,)2/3 * 



2/3^-2/3,4/3 



f^'Me, (25) 



corresponding to the mass ratio 

2/3 



Mrr, 



: 0.62 



77l/2(i?./i?t,)2/3 



. -1/3^-2/3 1/3 



(26) 



The core-to -envelope mass ratio increases with time, rapidly 
at first, then more gradually. The core radius, readily ob- 
tained from eq. (|20p . behaves in a qualitatively similar way. 



4 NUCLEAR BURNING 

4.1 "Main Sequence" Lifetime 

Fully convective stars are thought to be well-mixed, and 
therefore capable of burning a large fraction of their hydro- 
gen on the main sequence. Assuming an H mass fraction 
X = 0.75, the timescale to deplete all the hydrogen in a 
star of constant mass, radiating at the Eddington limit, is 
2 X lO'^ yr. If the star is growing at a steady rate by ac- 
cretion then the lifetime is doubled, to 4 x 10® yr. For a 
supermassive star growing by accretion, however, there are 
two factors that can shorten the lifetime if the star does not 
become fully convective. First, most of the mass may re- 
side in the convectively stable envelope, which does not mix 
with the core, thus restricting the fuel supply. Second, the 
thermonuclear luminosity required to maintain equilibrium 
exceeds the Eddington limit, because it must neutralize the 
binding energy of the freshly accreted matter in addition to 
replacing the lost radiation. The required luminosity is 



8 4nGMc 



-t 



(27) 



for a partially convective star with Afm/Af* > 0.1, and « 
Le for a fully convective star. To obtain the total nuclear 
energy required over the lifetime of the star, we integrate 
Lnuc over time. Comparing the total energy output with the 
fuel supply contained in the core, we obtain an estimate of 
the nuclear burning timescale, 

tnuc « 1.5 X 10'':^ yr. (28) 

We can now estimate the conditions inside a partially 
convective star at the end of hydrogen burning, and de- 
termine whether the star becomes fully convective. Sub- 
stituting eq. (dH) into eq. with x'^l^/rj^^^ « 0.9 and 



^^0.9m-/X8. (29) 

The core mass is independent of the accretion rate, 

« 1.2 X lO'^T-g^ Mq, (30) 
while the total mass is 

M, ^ 1.4 X lO^rhy^Tjyg Mq, (31) 

According to these estimates, a supermassive star will burn 
up its core hydrogen before reaching a fully convective state 
if m, > 0.7T-^. 



4.2 Central Temperature 

We have shown that the thermal relaxation of a growing 
supermassive star is sensitive to the central temperature of 
the convective core, which we have assumed to be constant 
in anticipation of the usual nuclear thermostatic effect. Up 
to now we have retained this temperature as a parameter: 
let us now estimate it for the dominant reactions likely to 
be occurring in these systems. 

The only thermonuclear energy source that can sup- 
port a supermassive star is the CNO cycle (Hoyle & Fowler 
1963a). At temperatures of interest the cycle reaches equi- 
librium quickly, hence we use the equilibrium reaction rate 
(Clayton 1983), assuming a hydrogen mass fraction X — 
0.75 but leaving the CNO mass fraction Zcno as a pa- 
rameter. We have also considered the possible role of the 
"hot" CNO mode (Mathews & Dietrich 1984) and the high- 
temperature saturation of the energy generation rate (e.g., 
Narayan & Heyl 2003), but have found both of these ef- 
fects to be unimportant at temperatures present in the core, 
due to the very low central densities of supermassive stars 
(p. = 0.04r,^/4m-|f T^g g cm"^). 

The energy generation rate within the convective core 
can be written in the form 



A 



™l/2y7/3 



n 



3/4 



6*3 (a;) exp 



B 



dx, (32) 



R*/Rt: 



0.6, we find 



where A and B axe constants. For a fully convective star, use 
rj = 4.8, replace m^.e by m.^g, and set Xm = 6.89685. In 
addition to computing the nuclear energy generation in the 
core, we will need to check whether any significant energy 
generation occurs in the envelope. 

To determine the equilibrium central temperature, we 
set Lc equal to the required energy generation rate given by 
eq. We use A = 1.10 x IO'^'^^cno erg s'^ B = 32.81 in 
eq. (|32|l to obtain I/c,cno- For all relevant values of Mm/M,, 
we can verify that essentially all nuclear burning is confined 
to the core. 

Figure 5 shows the equilibrium value of Tc,s as a func- 
tion of logifl(ZcNo/'Tiye ). For m*,6 ~ 0(1) and Zcno close 
to solar abundances (~ 0.01), we reproduce the result from 
Hoyle & Fowler (1963a) that T^g 0.7. Low tempera- 
tures favor more thermally relaxed stars (larger Mm/Mt) 
for a given accretion rate and elapsed time, primarily be- 
cause the core radius is inversely proportional to Tc- On the 
other, hand, for low CNO abundances (< lO^**) that may 
be present at high redshifts where supermassive stars might 
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Figure 5. Central temperature of the convective core in units of 

10** K, as a function of the CNO mass fraction Zq^q divided by 
1 /2 

g . Top, middle, and bottom curves are for envelope-to-core 
mass ratios Mf/Mm = 3, 10'^, and 10, respectively. (Note that 
Tc,s is not a monotonic function of M^/Mm-) 



be forming, Tc can exceed 10* K, possibly exceeding 2 x 10* 
under extreme conditions of low metallicity. Such high tem- 
peratures would imply that supermassive stars are never 
thermally relaxed, even for modest accretion rates < 1 Mq 

yr-^ 

If the star is totally devoid of CNO nuclei, hydrogen 
burns primarily via the PP III branch of the proton-proton 
chain (Clayton 1983). Accounting for neutrino losses, we 
have A = 3.71 x 10''^ erg s"\ B = 7.28 in eq. §2i- It is 
straightforward to show that the proton-proton chain can 
never meet the energy requirements of a supermassive star, 
nor can helium burning via the triple-a process. 

One can debate the likelihood that a supermassive star 
with such low CNO abundance (^cno ^ 10^*) could ever 
form under cosmological conditions (e.g., Trenti & Stiavelli 
2009; Trenti, Stiavelh & ShuU 2009). The point may be 
moot, however, because a supermassive star will create its 
own metals. To see this, we estimate the production rate of 
CNO elements via the triple-a process in the core: 

McNO = 



X 6g (x) exp 



44.03 

'n,803{x) 



dxMcc yr"V(33) 



At Tc,8 ~ 2, the core produces ~ 1O~^M0 yr~^, and at 
Tc,8 ~ 3 the CNO production rate is increased > 1O~^M0 
yr~^. If the star can persist for a few thousand years in 
Kelvin contraction at Tc.s < 3, it should be able to produce 
enough CNO elements to support itself by burning hydro- 
gen. 



5 BLACK HOLE FORMATION AND GROWTH 

5.1 Formation of a Seed Black Hole 

After exhausting its hydrogen, the core of a supermassive 
star will contract and heat up until it suffers catastrophic 
neutrino losses and collapses. If the mass shell M encloses 
a specific angular momentum less than ~ GM/c, it can col- 
lapse directly to a black hole without any angular momen- 
tum transport. Approximating the convective core as a uni- 
form sphere of density pc and radius Rm, we can estimate 
the maximum ratio of angular velocity, Q., to local Keplerian 
angular velocity, Q.k = {GM / R^Y^^ , consistent with direct 
black hole formation: 



< 



GM{R) 
c^R 



Qmm]ilT.,,(^-^) . (34) 



Note that this constraint becomes more stringent at smaller 
radii. If, for example, the convection drives the core to solid 
body rotation, Q, ~ const., then e{R) is also constant and 
the core might satisfy condition (|34p at large radii but not 
near the center. Given that the supermassive star is formed 
from gas that probably has substantial rotation, this strin- 
gent constraint suggests that the formation of the black hole 
depends on angular momentum transport in the collapsing 
gas. Since the collapsing gas within any mass shell is strongly 
self-gravitating, angular momentum transport can be driven 
efficiently by global gravitational torques and turbulence re- 
sulting from nonaxisymmetric gravitational instabilities, as 
discussed in Begelman et al. (2006) — although magnetic 
torques could also be important. 

The transport of angular momentum outward is un- 
avoidably accompanied by the outward transport of energy 
(Blandford & Begelman 1999) . An unknown fraction of this 
liberated energy might escape the system via a jet punching 
through the star (as in the collapsar model of gamma-ray 
bursts). Neutrino losses might also be important, particu- 
larly when the black hole mass is small. The remainder of the 
energy must pass through the star on its way out, presum- 
ably driving a strong convective flux through the outer core 
and surrounding envelope (Blandford & Begelman 2004). In 
associating the convective luminosity with the rate of black 
hole growth, Lbh = cMbhc^, we fold in these possible in- 
efficiencies through the factor e. We recognize the standard 
efficiency of black hole accretion by using the normalization 
e = 0.1e_i, with the understanding that e_i might be ^ 1. 

Initially, the liberated energy flux is trapped in the 
star, causing it to inflate. As we show below, this hap- 
pens so quickly that h'h can be taken to be constant dur- 
ing this process. Since R oc M^^'^T~^ for a radiation- 
dominated, convective body, the core temperature drops 
rapidly as the star expands, as does the total binding en- 
ergy. An upper limit to the energy that can be absorbed 
without dispersing the star is the initial binding energy, 
Eb = Eg -\- Erot (neglecting general relativistic corrections) . 
For a fully convective star 
while Erot = 3.4 x 10' //i, g 
the "potential energy-weighted" mean rotation parameter 
J eGM{R)dM/R. This leads to a maximum "seed" 
black hole mass that can grow during the intial expansion 
phase: 



2.3 X 10''='m.,6rc,8 erg 
^''ml^RTc,s{e)w erg, where {e)w is 



Meocd = = 130 



(^l + 150ml(^{£)w^ Mq. (35) 
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The estimate of Ef, (and thus Mgeed) is easily generalized to 
the case where the supermassive star is not fully convective 
when the black hole starts to grow. 

To estimate the timescale for this initial growth, we 
apply the modified Bondi accretion model described in Be- 
galman et al. (2006, 2008). Because the accretion radius of 
the seed black hole is small compared to the radius of the 
star, we can use the central pressure and density as approx- 
imate outer boundary conditions for the flow. The standard 
Bondi (1952) formula for the accretion rate is reduced by 
a factor ~ e~^{ca/c}'^, where Cs = (4Pc/3pc)^''^ is the adi- 
abatic sound speed outside the accretion radius. This fac- 
tor is necessary to ensure that the convective flux escap- 
ing from the accretion flow does not exceed the saturated 
value ~ PcCc (Gruzinov 1998; Blandford & Begelman 1999; 
Narayan, Igumenshchev & Abramowicz 2000; Quataert & 
Gruzinov 2000). For a black hole mass of mBnAf©, the 
growth time is 



where rriifi ~ min 



, MbH -1 3/4^^-5/2 

ten = -T— = 450 e-lmBH"^.,6 T^.g yr, 



(36) 



implying that a black hole of modest initial mass will quickly 

grow to Msecd. 

The growth of the black hole cannot completely disrupt 
the star, because this would interrupt the fuel supply and 
stop accretion. Once the black hole reaches a mass Msecd, 
the convective envelope expands at approximately constant 
mass, in response to the black hole power 

r r, 1^42 2 -3/4rrn5/2 -1 D-5/2 /oryS 

Lbh = 2 X 10 e-imseod"i, 6 T^/g erg s cc ' . (37) 

The star continues to expand as long as Lbh exceeds the 
Eddington limit for the star, Le = 1.4 x lO^^Mt^e erg s"^. 
Once the limit is reach, the black hole-star system should 
come into a kind of equilibrium that we have dubbed a "qua- 
sistar" — a bloated convective envelope powered by black 
hole accretion at the center. 



1/2^ 



is the initial mass of 



4m«,1.4m, j.^^ 

the quasistar in units of 10^ M©, from eq. (I3ip . 

From eq. (11) of Begelman et al. (2008) (with parame- 
ters set to their flducial values), the effective temperature of 
the quasistar's photosphere (in units of 10'^ K) is given by 



Tph,3 = 2 X 10 e: 



-1/5 -2/5 7/20 



(39) 



which decreases with time as the black hole and quasis- 
tar grow (Begelman et al. 2008). The quasistar disperses 
and black hole growth stops when Tph reaches Tmin ~ 4000 
K. In contrast to a contracting protostar or red giant en- 
velope which manages to maintain dynamical equilibrium 
while evolving at close to the minimum effective tempera- 
ture (the "Hayashi track"), the quasistar loses equilibrium 
and disperses at this point because of the nature of the black 
hole energy source as well as the fact that it is supported by 
radiation pressure rather than gas pressure. This behavior 
is discussed at length in Begelman et al. (2008). 

Given the estimates above for Afsecd and Mi, it is clear 
that the black hole is generally able grow by a large factor 
before the quasistar disperses. Assuming that the final black 
hole mass, ?n/_BH is much larger than rriscad, we find 



m/,BH 



1.3 X 10^ /4600 K\^/2 

f ' \ -t- mm / 



2 X 10^ 7/8 

m , 



1/2 



f,6 ■ 



(40) 



where Mf is the final mass of the quasistar. For modest 
accretion rates m« < 0(1), the quasistar does not grow 
much between its formation and dispersal, Mf « Mi, and 
for m, ~ 0{1) the black hole left behind when the quasistar 
disperses has a mass in the range ~ (10* - 105)£:J/^Mq. 
But for m, ^ Se^^T^'^g^ , the final quasistar mass is consid- 
erably larger, m/_6 ~ lAet/imV^ and the final black hole 

o w in4 -1/9 ■ 7/9 

mass IS m/,BH 3 x 10 e_j m* . 



5.2 Growth Inside Quasistar 

As discussed in previous publications (Begelman et al. 2006, 
2008), under suitable initial conditions a black hole can grow 
rapidly and by a large factor inside a quasistar. The growth 
occurs at a highly super-Eddington rate (i.e., with mass e- 
folding in much less than the Salpeter time, ~ 40 Myr) be- 
cause the accretion power of the black hole is regulated at 
the Eddington limit corresponding to M, , which is much 
larger than Mbh. The main requirement for stable growth is 
that the ratio of black- hole to envelope mass be smaller than 
~ 0.01 (with some dependence on envelope abundances, via 
the opacity); otherwise the effective temperature becomes 
so low that the envelope cannot radiate away the flux, and 
it disperses. We refer readers to Begelman et al. (2008) for 
detailed discussion of quasistar structure and evolution. 

To calculate the co-evolution of the black hole and qua- 
sistar, we apply the model described in section 4 of Begel- 
man et al. (2008), with the two feedback efficiency parame- 
ters, a and e, set equal for simplicity. Integrating eq. (51) of 
Begelman et al. (2008) and using the quasistar mass as the 
independent variable rather than time, we obtain 

1.3 X 10* .2 2 s 

mBH = Jriaccd H : (m, g-mig), (38) 

e_im, ^ ' 



6 SUMMARY AND DISCUSSION 

We have revisited the theory of supermassive stars, tak- 
ing into account the fact that their short nuclear burning 
timescales {^4 Myr) imply that the matter forming them 
must accumulate rapidly. If the accumulation of matter is 
rapid enough (typically, more than a few solar masses per 
year), then the star never has time to become thermally re- 
laxed and fully convective. Instead, its structure consists of 
a convective core containing a minority of the mass and oc- 
cupying a small fraction of the stellar volume, surrounded by 
a convectively stable envelope with the specific entropy in- 
creasing outward as enclosed mass to the 2/3 power. We have 
calculated the structures of these envelopes — for which we 
propose the term "hylotropes" (following a suggestion by 
A. Accardi and G. Lodato) — and the conditions for match- 
ing them to the convective cores, and have thus devised a 
theory for the evolution of partially convective supermassive 
stars. 

Partially convective supermassive stars are strikingly 
different from their fully convective counterparts (described 
by n = 3 polytropes) in a number of respects. They do 
not follow the same mass-radius relation as fully convective 
stars. Instead of the scaling R oc M^^'^, which applies for 
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a fully convective supermassive star with a central temper- 
ature regulated by thermonuclear reactions, our partially 
convective models have a radius proportional to the (as- 
sumed constant) accretion rate. Physically, the stellar radius 
is comparable to the "trapping radius" within which radia- 
tion is unable to escape inward advection by the accumulat- 
ing gas (Begelman 1978, 1979). Also, their nuclear-burning 
("main sequence") lifetimes can be shorter because the fuel 
in the convectively stable envelope is never mixed into the 
core. 

Hydrogen-burning supermassive stars in the mass range 
considered here (> lO^M©) can exist only if stabilized by ro- 
tation or some other "stiff' form of energy (such as magnetic 
fields or turbulence) (Fowler 1964, 1966). The minimal level 
of rotation required to cancel the general-relativistic insta- 
bility is not dynamically significant. However, the angular 
momentum may well be high enough to inhibit the forma- 
tion of the black hole, once the core runs out of hydrogen 
and collapses. We have argued that the initial mass of the 
seed black hole, if limited by the energy released in forming 
it, may be as small as a few hundred solar masses — i.e., 
only 10"* of the total stellar mass. However, there is likely 
to be a second phase of rapid black hole growth once the 
inflated stellar envelope reaches equilibrium as a red giant- 
like quasistar. We therefore estimate that supermassive stars 
should leave behind black holes with typical masses of a few 
percent that of the star, i.e., roug hly ~ 10" - 10^ Mq. 

These estimates are admittedly quite uncertain, mainly 
because we do not understand how energy liberated by black 
hole accretion couples to the remainder of the star. We can- 
not rule out the possibility that most of the energy punches 
through the star in a pair of jets (in a fashion similar to 
the coUapsar model of gamma-ray bursts [Woosley 1993]), 
in which case both the seed and final black hole masses could 
be much larger. We note that while our estimate of the ini- 
tial black hole mass is qualitatively similar to values (few 
tens of Mq) predicted by Begelman et al. (2006), our new 
value for the predicted envelope mass at the time of black 
hole formation is 2-3 orders of magnitude larger. The rea- 
son for this discrepancy is that Begelman et al. (2006) failed 
to account for the rapid growth of the convective core due 
to heating by efficient thermonuclear reactions in the CNO 
cycle. (Effectively, we considered pure Pop HI abundances, 
for which thermonuclear energy release is too weak to stall 
the formation of the seed hole.) 

Supermassive stars can be strong sources of far-UV 
radiation, provided that their spectra are not severely de- 
graded by reprocessing at radii beyond Rti . Taking the pho- 
tospheric radius to be 0.67?tr, we find an effective tempera- 
ture TcH = 1.2 X lO^ml^^rh^^^^ K for a partially convective 

star of mass M., corresponding to 1.3 x 10^ rh'^^^T'^^'^ K 
at the end of hydrogen-burning. However, the infalling gas 
joining a supermassive star is already very optically thick at 
Rtr (t ~ c/v, where v is the infall speed), implying that most 
of the radiation escapes from further out and is likely to be 
softer. On the other hand, the opacity is strongly scattering- 
dominated, so the color temperature could be several times 
higher than the effective temperature. The hardest radiation 
is likely to escape from the polar regions, where rotational 
effects decrease the density. 

Assuming a blackbody spectrum with the effective tem- 



perature estimated above and m, = Tc.s ~ 1, we find 
that most of the radiation is capable of ionizing hydrogen 
(E > 13.6 eV), with an ionizing luminosity Lion ~ 2 x 10** 
erg s~^ and photon flux A/ion ~ 4 x 10^* s^^. The star ra- 
diates into the Lyman-Werner band (11.2 < E < 13.6 eV), 
capable of dissociating H2, with Ldiss ~ 5 x lO*'^ erg s~^ 
and Adiss ~ 3 x 10^"^ s~^. To produce an interesting mean 
intensity of dissociating radiation, J„ > 10~^^ erg cm~^ s~^ 
Hz^^ sr~^ between z ~ 12 and z ~ 6, would require (very 
roughly) a comoving density of supermassive stars of at least 
~ 10~* Mpc~^. Given the short lifetimes of these objects, 
it is not clear how plausible this is. We will investigate cos- 
mological scenarios for the formation of supermassive stars 
in a future publication. 
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